Anisotropic Transport Properties of Ferromagnetic-Superconducting Bilayers. 
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We study the transport properties of vortex matter in a superconducting thin film separated by a 
thin insulator layer from a ferromagnetic layer. We assume an alternating stripe structure for both 
FM and SC layers as found in 0. We calculate the periodic pinning force in the stripe structure 
resulting from a highly inhomogeneous distribution of the vortices and antivortices. We show that 
the transport in SC-FM bilayer is highly anisotropic. In the absence of random pinning it displays 
a finite resistance for the current perpendicular to stripe and is superconducting for the current 
parallel to stripes. The average vortex velocity, electric field due to the vortex motion, Josephson 
frequency and higher harmonics of the vortex oscillatory motion are calculated. 
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The interest in Heterogeneous ferromagnetic- 
superconducting systems has grown rapidly in recent 
years. This interest stems not only from their possible 
technological applications but also from new physical 
phenomena arising from the interaction between two 
order parameters. Typically such a system consists of 
a superconductor (SC) placed in close proximity with a 
periodic ferromagnetic structure (FS) such as an array 
of ferromagnetic dots or holes. The two systems are sep- 
arated by an infinitely thin layer of insulator oxide that 
guarantees the suppression of proximity effects. It was 
demonestrated experimentally [lj-|3j that the interaction 
between the superconductor and the ferromagnet may 
lead to formation of superconducting vortices interacting 
with the FM film. Theoretical studies of such systems 
have been done in |j]-0 

Recently, Erdin ei.aL0 studied the equilibrium struc- 
ture of a FM-SC bilayer (FSB). They have proved that 
it represents a two-dimensional periodic stripe domain 
structure consisting of two equivalent sub-lattices, in 
which both the magnetization m z (r) and the vortex den- 
sity n v (r) alternate. Thus, they predicted spontaneous 
violation of the translational and rotational symmetry in 
the bilayer. In this article we study the transport prop- 
erties of the FSB. They are associated with the driv- 
ing force acting on the vortex lattice from an external 
electric current. We show that the FSB exhibits strong 
anisotropy of the transport properties: the bilayer may 
be superconducting for the current parallel to the domain 
walls and resistive when the current is perpendicular to 
them. 

Periodic pinning forces in the direction parallel to the 
stripes do not appear in continuously distributed vor- 
tices, their reappearance is associated with the discrete- 
ness of the vortex lattice. Therefore, we need to modify 
the theory to incorporate the discreteness effects. Let 
us assume that the saturation magnetization per unit 
area of the FM film is m and its width is L. The en- 
ergy necessary to create a single Pearl vortex || in the 



superconductor is e v o = eoln(^), with eo = 16 y?j A where 
(j>a is the flux quantum, A 
tion depth 9], Xl 



is the effective penetra- 



4 ~ C e 1 ^ s ^ ne London penetration 
depth, d s is the thickness of the superconducting layer 
and £ is its coherence length. It was shown in [?J that the 
interaction between the superconducting vortices and the 
magnetization in the stripe structure renormalizes the 
single-vortex energy to the value e v = e v0 — m0 o which 
must be negative to allow development of the stripes. 
The density of the superconducting vortices increases 
when approaching the domain walls and can be expressed 



as n v (x) 



where fa 



is the renor- 



L0 O sin(^) 

malized magnetization of the FM stripe. The vortices 
spontaneously appear in the superconductor. We assume 
that the vortices inside one stripe are arranged in paral- 
lel chains. Each chain is periodic with the same lattice 
constant b along the chain, whereas the distance between 
fc-th and (fc+l)-th chain ak depends on k. The correspon- 
dence between this discrete arrangement and continuous 
approximation is established by the requirement that 
the local vortex density n v (xk) calculated in must be 
equal to (&a&) _1 . The coordinate Xk is determined in 

terms of ak as a sum: xu — 2fc'=o flfe ' ■ For definiteness 
we choose the origin in the center of the stripe. We as- 
sume that the total number of the vortex chains 27V in 
a stripe is large. Then some of them are located very 
close to the domain walls. Let us remind that, in the 
continuous approximation n v = ^rm\ ; where L 
is the domain width. Considering the nearest to the do- 
main wall vortex chain (with the number, N), we put 
n v {x?f) = On the other hand x^ = L — a^ . Since 
2^ < 1, we find: b = The total number of chains 

L ' m 

in a stripe is 27V, where N = b f Q n v (x)dx = | ln(j)- 
We cut the integration (and summation) at a distance 
~ A from the domain wall where the continuous approxi- 
mation breaks. Thus the minimum value of a is A. When 
transport current passes through the superconducting 
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FIG. 1: Schematical vortex distribution in the FM-SC bilayer. 
The sign ± refers to the vorticity of the trapped flux. 



film, the vortices start to move. To simplify the problem 
we assume that all vortices in each stripe move together 
as well as all antivortices in the neighboring stripe do. We 
denote their positions r + = (x + ,y + ) and r_ = (a;_,y_), 
respectively. Forces acting on a moving vortex are the 
Magnus force, the viscous force and the periodic pinning 
force. The Magnus force is f m = Trn s fid s (v s — r) x z, where 
n s is the superconducting electron density, v s is the ve- 
locity of the superconducting electron and r is the vortex 
velocity. The viscous (friction) force is f/ = — rjr where 



PnC* 



is the Bardeen-Stephen drag coefficient [H 



V 

H C 2 is the upper critical magnetic field, p n is the resis- 
tivity of the superconducting sample in the normal state, 
and c is the speed of light. The periodic pinning forces 
are due to the interaction of the vortex with the pinning 
centers and the domain walls. In the FM-SC bilayer the 
pinning force is due to the interaction of the domain walls 
with the vortices and antivortices and the vortex-vortex 
interaction U vv given by: 



1 



rc„(r)F(r-r')n„(r')<2 2 rcfV 



(1) 



wher V(r — r') is the pair interaction between vortex 
located at r and another at r'. When |r — r'| 3> A the 

pair interaction can be written as V(r — r') = j^- rrzpi 
The interaction energy between two parallel chains 
located at xi and xy and vertically shifted with respect 
to each other by an interval bs (s <C 0) reads: 



U(xi,xi',s) 



N 

E 

n,m— 



1 y/ (xi - XI') 2 + (n - 



2 b 2 



(2) 



where Nq is the number of vortices or antivortices in a 
single chain. For infinite chains (No — > oo) Eq.J2J) can be 
rewritten as 



U(x h x v ,s) 



E 

k= — oo 



No4> 2 „ 



y/(xi-Xl') 2 + (k- S ) 2 b 2 



(3) 



The sum in can be calculated using Poisson summa- 
tion formula |l2j| . Since the force is zero in the continuous 
approximation, it is possible to retain the lowest non-zero 
harmonic in the Poisson summation. Thus, we arrive at 
the following interaction energy of two chains: 



U(x l ,x v ,s) = 



4n 2 b 



cos(2tts)xh' 



(4) 



where \u 
chains la;; 



e "" b . The distance between two 
xi'\ is larger or equals A hence xiv ~ X — 
e _J f-. Typical value of xw is e~^~ where S m = — 
gS n ™ d ^ , with g is Lande factor, S is the ferromagnet el- 
ementary spin, n m and d m are the electrons density and 
thickness of the magnetic film respectively. 

We conclude that the amplitude of the periodic po- 
tential for displacements parallel to the domains is ex- 
ponentially small in units of eo, near the transition tem- 
perature. Relative displacements in perpendicular direc- 
tion have energy barrier ~ eo even in continuous ap- 
proximation. We model the restoring forces by sim- 
ple sines dependencies f x = f±sin( 2 ^-(x + — X-)), f y = 
/I, sm(^(y+ - !,_)), where /j_ ~ a an d /„ ~ fe'^ « 
fx- 

When the supercurrent is perpendicular to domains, 
the equations of motion for a vortex and antivortex 



T)X + 



F 

F x+ 



/||sm( — (y+ 



y-)) 



—y+ + fx sin — (x+ - x. 

v s a 



)) 



vv- = 



r/X- = 



F 

-F+—X+ 



f\\ sin (y(yn 



-y-)) 



V- - fx sm — (x+ - X-)) 

v. a 



where F = nn s hd s v s , fx 



(5) 

(6) 

(7) 

(8) 
If the 



IX = f and /|| = 
current is smaller than a critical value, Eq. l">|St accept a 
static solution 



X-\- 



y+ 



Fb 
b 



-y- 



arcsm( — j 
f\\ 

arcsm(— 
4tt /|| 



(9) 
(10) 



It is valid at F < f\\ = ^f 2 -. For F > fit or equivalently, 
if the current is larger than its critical value, vortices and 
antivortices start to move. The solution of Eg .1 151811 for 

F > F c reads: 



Xj^- X_ 
X + + X — 

y+ - v- 
y+ + v- 



= o 



—(y+-y- 

r)V s 

^arctan^ 




(11) 
(12) 



1- JLtan(u;£t))(13) 
(14) 
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The time-averaged local field components are 



FIG. 2: The vortex displacement as a function of time in the 
overcritical regime. Time is measured in units of — K- and j/+ 

in units of b and x = 10 -4 . 



where ujq 



— bi( Jp-i +v i v ij is the Josephson frequency. 



Thus the vortices and antivortices acquire the same ve- 
locity components v 



+x 



<j- x in the direction of the 
current and opposite velocity components v +y = —V- y 
in the direction perpendicular to the current. The do- 
main walls do not interfere such a motion if they move 
in the direction of the current with the same velocity 
Vdw = v+x = v~x as vortices and antivortices. Such a mo- 
tion is a Goldstone mode. The solution (|1 1113(1 displays 
an oscillatory motion of the vortices and antivortices in 
the direction parallel to the domain walls, in addition to 
their motion together with the domain walls along the 
direction of the current. Higher harmonics of the vortex 
motion can be calculated analytically. The distribution 
of vortices (antivortices) is inhomogeneous in the direc- 
tion perpendicular to the domains. The local electric field 
E due to the vortex motion is related to its time-average 
velocity < v+ > as E = -^-n v (r){< v + > xz) [Hf. 
Therefore, the local field produced by vortices in the di- 
rection parallel to the domains is equal but opposite in 
sign to the one produced by antivortices, while the local 
field produced by vortices and antivortices in the direc- 
tion perpendicular to the domains has both equal magni- 
tude and sign. The time average of the vortex (antivor- 
tex) velocity over a period T = ■^r is 



= ±- 



(V 2 + W) 



< v +x > 



F < V+y > 

rjv s 



(15) 
(16) 



E , 



Ey=T 



mF 



F 2 -L 



(17) 
(18) 



The upper sign in Ea. l(15|) and Ea. l(18|) refer to the vor- 
tices velocity and produced field along the domain while 
the lower sign refer for those due to antivortices. Non- 
zero average electric field due to all vortices and antivor- 
tices in the FSB appears only in the direction perpendicu- 
lar to the domains. The critical current J c is related to F c 
as J c = -T^j-F c . Plugging F c — £ ° < ^ x into the expression 
for J c and accepting x = 10~ 4 -10~ 2 , b = 10~ 4 -10~ 5 cto, 



and 



10 22 cm" 3 we find: J r ~ 10 3 - 10 



5 A 



When the current flows parallel to the stripes, the 
FM domain walls stay at rest while vortices and an- 
tivortices move both parallel and perpendicular to the 
domains. The solution of equations of motion for vor- 
tices and antivortices shows that they move opposite 
to one another both in x and y directions. Their mo- 
tion along x is oscillatory with fundamental frequency 

ujj = 2 ^HF' i W : ^Y 2 '- Tne mot i on of vortices and an- 
tivortices in the parallel direction proceeds until the dis- 
tance between them becomes half lattice spacing |. Once 
the vertical shift between the vortices and antivortices 

The critical current in 



reaches |, their motion freezes, 
this case is J c = , the lattice spacing a is of the 
— 10 -4 cm, hence the critical current 
10 7 - lO 8 ^, which is at least 10 2 



order of A ~ 1CT 
J c is of the order 
times larger than the critical current for parallel current. 
Therefore, the system may be superconducting for the 
current parallel to the stripes and exhibit finite resistance 
for perpendicular current. The difference in the critical 
currents for parallel and perpendicular directions is due 
to the exponential factor \ which is small if b <C A. The 
anisotropy is pronounced when S m is large which can be 
achieved by using thicker FM layers and decreasing the 
density of the superconducting electrons. 5 m is temper- 
ature dependent and eventually decreases when temper- 
ature decreases starting from T s . However, at the tem- 
perature of vortex disappearance T v < T s the value rh 
turns into zero and \ again becomes exponentially small. 
Thus anisotropy has a minimum between T v and T s . 

Kopnin and Vinokur considered a collection of su- 
perconducting grains with the washboard pinning poten- 
tials a model of random pinning. They obtained a simi- 
lar result for vortex sliding in external magnetic field with 
a supercurrent applied. In contrast to their work (they 
considered vortices only), we consider vortices and an- 
tivortices in the field of periodic pinning and completely 
neglect the random pinning. 

Let us discuss briefly how the magnetic field generated 
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by supercurrent changes our result. In [16j it was shown 
that at sufficiently small critical magnetic field the do- 
mains vanish. Therefore, in general, magnetic field sup- 
presses both the anisotropy and periodic pinning at a 
critical field for which domains disappear. At such criti- 
cal field only random pinning prevails. However, the total 
per unit length current is proportional to the thickness 
of the SC film and can be kept small. 

In conclusion, we studied the transport properties of 
the FM-SC bilayer in a state with stripe domains of al- 
ternating magnetization and vorticity. We showed that, 
in the absence of a driving force, the vortices and antivor- 
tices lines themselves up in straight chains configuration. 
We showed that the force between two chains of vortices 
falls off exponentially as a function of the distance sepa- 
rating the chains. We argued that the exponential decay 
of the pinning force in the direction parallel to the do- 
mains drops faster in the vicinity of the superconducting 
transition temperature T s and vortex disappearence tem- 
perature T v . We solved the equations of motion for vor- 
tices and antivortices in two cases for the driving current 
direction, parallel or perpendicular to the domains. We 
calculated the critical current for both cases and found 
that its value for the parallel current is much higher than 
for perpendicular one. Our most important result is a 
strong anisotropy of the critical current. We expect the 
ratio of the parallel to perpendicular critical current is in 
the range 10 2 -j- 10 4 close to the superconducting tran- 
sition temperature T s and to the vortex disappearance 
temperature T v . The anisotropy decreases rapidly when 
the temperature goes from the ends of this interval reach- 
ing its minimum somewhere inside it. This anisotropic 
transport behavior could serve as a diagnostic tool to 
discover spontaneous topological structures in magnetic- 
superconducting systems. 
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